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Abstract. In inhomogeneous optically thick synchrotron sources a substantial part of the electron population 
at low energies can be hidden by self-absorption and overpowered by high energy electrons in optically thin 
emission. These invisible electrons produce Faraday rotation and conversion, leaving their footprints in the linear 
and circular polarized radiation of the source. An important factor is also the magnetic field structure, which can 
be characterized in most cases by a global magnetic field and a turbulent component. We present the basic radiative 
transfer coefficients for polarized synchrotron radiation and apply them to the standard jet model for relativistic 
radio jets. The model can successfully explain the unusual circular and linear polarization of the Galactic Centre 
radio source Sgr A* and its sibling M81*. It also can account for the circular polarization found in jets of more 
luminous quasars and X-ray binaries. The high ratio of circular to linear polarization requires the presence of 
a significant fraction of hidden matter and low-energy electrons in these jets. The stable handedness of circular 
polarization requires stable global magnetic field components with non-vanishing magnetic flux along the jet, while 
the low degree of total polarization implies also a significant turbulent held. The most favoured magnetic held 
conhguration is that of a helix, while a purely toroidal held is unable to produce significant circular polarization. 

If connected to the magnetosphere of the black hole, the circular polarization and the jet direction determine 
the magnetic poles of the system which is stable over long periods of time. This may also have implications for 
possible magnetic held conhgurations in accretion hows. 

Key words. Polarization - Radiation mechanisms: non-thermal - Radiative transfer - Galaxies: jets - Accretion, 
accretion disks - Radio continuum: general 


1. Introduction 

The detection of circular polarization (CP) in the contin¬ 
uum of radio sources is believed to be a powerful tool to 
test physical source models (Hodge & Aller 197S| ) . But CP 
in extragalactic radio sources is extremely elusive (Roberts 
et al. 1975; Ryle & Brodie 1981; Weiler & de Pater 1983) 
and is detected in only a few sources. 

A more recent ATCA-survey (Rayner, Norris & Sault 
2000| ) for CP in radio-loud Quasars, BL Lacs and Radio 
Galaxies with improved sensitivity of 0.01%, has shown a 
clear correlation of fractional CP with spectral index, in 
the sense that CP is stronger in flat and inverted spectrum 
sources. Circularly polarized radiation is therefore prefer¬ 
entially produced in self-absorbed radio cores. The frac¬ 
tional CP at 5 GHz is found to be between 0.05% and 0.5% 
in 11 out of 13 inverted spectrum sources at the ATCA 
spatial resolution of 2 arcsec. At hig her VLBA-resolution 
(~ 0.5mas) Homan & Wardle ( 1999|) report localized CP 
of 0.3%-l% in the jet-cores of 3C273, PKS 0528+134, and 
3C279, which in a few cases may be as high as the local 
linear polarization. It is also found, that intraday variable 
sources are circularly polarized (Macquart et al. 2000|) , 


and that LP (linear polarization) and CP are both variable 
on timescales below 1 day. Recently CP was also found in 
X-ray binaries (Fender et al. 2000| fc2002). 


While the handedness of CP is remarkably stable over 


several years (Komesaroff et al. 1984; Homan & Wardle 
|1999 ; Fender et al. [2002| ) for individual sources, no pre¬ 
ferred handedness of CP in general is found. 

An even more challenging situation than observed in 
radio-loud extragalactic jet sources presents itself in the 
centre of our galaxy. The compact radio source Sgr A* 
(see Melia & Falcke [2001 ), believed to be coincident with 
the dynamical centre of the Milky Way with a mass of 
2.610 6 Mq (Eckart & Genzel 199G; Ghez et al. 1998) pre¬ 
sumably in a single black hole, exhibits consistently larger 
circular than linear polarization in the range of 1.4 to 15 
GHz (Bower et al. 1999b; Sault & Macquart 1999| ) with 
CP between 0.2% and 1%. LP is small and below the de¬ 


tection limits (Bower et al. 1999a , 1999c ) with the excep¬ 
tion of sub-mm measurements, which possibly shows LP 
at a l evel of 10% in the range 750 — 2000/rm (Aitken et al. 
|2000 ). The beam size of the sub-mm observations is ~ 10 
arcsec. The flux is dominated by extended dust emission 
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or free-free emission and the synchrotron source is com¬ 
parably weak at these wavelength. 

The inverted radio spectrum of Sgr A* (S u oc v a , a ~ 
0.3) can be interpreted as either optically thin synchrotron 
emission (Beckert et al. 1996) or self-absorbed synchrotron 
emission from a jet-like outflow (Falcke et al. 1993; Falcke 
& Markoff [2000 ). The idea of synchrotron emission by 
thermal electrons from Sgr A* was briefly considered by 
Reynolds & McKee (198C) and revived for mildly rela¬ 
tivistic electrons in the self-absorbed ADAF models for 
accretion in the galactic centre (Narayan et al. 1998). The 
first ADAF models under-predicted the radio flux between 
1 100 GHz, which can be attributed to an outflow or jet. 
The upper limits for Sgr A* in the infrared require a sharp 
high energy cut-off for the electron distribution below 7 max 
of a few xlO 2 . Therefore thermal or quasi-monoenergetic 
electrons are re sponsible for the radio emission (Beckert 
& Duschl 1997 ), which distinguishes Sgr A* from high- 
luminosity, radio-loud AGNs. A close relative of Sgr A* 
is found in the centre of the normal spiral M81. The ra¬ 
dio source M81* e xhibit s an elongated jet-like structure 
(Bietenholz et al. 2000), has a similiar radio spectrum 
(Reuter & Lesch, 1996 ), a slightly larger luminosity, still 
below the AGN level, and has recently be found to be 
circularly polarized (Brunthaler et al. [2001 ) without de¬ 
tectable LP. 

The fractional variability of CP is usually stronger 
than of LP, which in turn is stronger than for the to¬ 
tal intensity. Together with the preserved handedness this 
poses servere constrains on possible scenarios for CP pro¬ 
duction and its variability (Komesaroff et al. 1984). The 
suggested mechanisms are (a) intrinsic cyclo-synchrotron 
emission from low-energy electrons or from electrons with 
small pitch angles see n clos e to the magnetic field direc¬ 
tion (Legg & Westfold 1968 ), conversion from LP to CP as 
a propagation effect induced by (&) low ener gy el ectrons 
inside the relativistic plasma (Hodge & Aller 1977 ) or (c) 
by a magnetized cold plasma surrounding the synchrotron 
source. This requires either Faraday rotation (not possi¬ 
ble in pure electron/positron jets) or changing (e.g., turbu¬ 
lent) B-field directions along the line of sight in the source. 
A further possiblity for CP production are ( d ) inhomoge¬ 
neous rotation measures in intervening cold plasma either 
close to the source or in our galaxy (Macquart & Melrose 
200C| ) . The existence of these plasma screens can be infered 
from interstellar scattering believed to be the cause for 
intraday variability in some sources (Rickett et al. |1995; 


Dennett-Thorpe & de Bruyn 2000; Macquart et al. 2000 


Beckert et al. 2002). This model predicts variable CP with 
a time averaged mean of <CP>= 0. 

In this paper we consider propagation effects like 
Faraday rotation and cyclic conversion of LP to CP and 
back (Pacholczyk 1973| ) in turbulent, self-absorbed jets or 
outf lows. First results were already published in Falcke et 
al. ( 2002| ). We rederive some of the basic radiation trans¬ 
fer coefficients which, for example, could also be used for 
anisotropic particle distributions. The application of con¬ 
version to compact radio jets has been explored perviously 


by Jones ( |1988 ) using different techniques and without 
focusing on sources with large circular polarization and 
the role of globally ordered magnetic fields. Here we in¬ 
vestigate the standard jet model with respect to the new 
polarization data placing some emphasis on the role of 
turbulence, the ratio of low- to high-energy particles, and 
the magnetic field confirguation. 

The paper is organized as follows: In Sec. |] we re¬ 
view the basic production channels for CP. The outfow/jet 
model and the possible turbulence in the 13-field is pre¬ 
sented in Sec.|. The consequences of Faraday rotation 
and conversion are discussed in Sec.f| followed by a de¬ 
tailed model of Sgr A*. Polarization variability is the topic 
of Sec. H and we close with a discussion of our results in 
Sec. |[ 


2. Polarized Synchrotron Emission & Radiative 
Transfer 

2.1. Synchrotron Emission 

Relativistic electrons or positrons gyrate in a magnetic 
field B with a frequency v gT = vb^Ii- The basic cyclotron 
frequency 


r'B± = vb sin if) 


qB sin ijj 
27rm e c 


= 2.8 • 10 6 Hz 


B 

1G 


sin^ 


( 1 ) 


depends on the field B± = B sin ip perpendicular to the 
particle momentum. The transition from cyclotron to syn¬ 
chrotron radiation occurs, when emission at higher mul¬ 
tiples of the gyro-frequency, which are weak for non- 
relativistic particles, become stronger than emission at u gI 
due to relativistic boosting in the instantaneous direction 
of the particle. The synchrotron spectrum of individual 
particles reaches its maximum at v c = l/2y 2 VB- A popu¬ 
lation of particles with distributions in energy and pitch 
angle ip produce an emission spectrum, in which the indi¬ 
vidual cyclotron lines are blended together and a smooth 
spectrum emerges. The cyclotron emission at low frequen¬ 
cies is circularly polarized and reflects the spiral motion 
of the particles. The emission at the maximum at v c is 
seen only for a time y -2 of the gyration period and the 
weakly curved motion during that interval produces the 
large linear polarization of synchrotron emission perpen¬ 
dicular to the B-field. The emission of individual electrons 
is also circularly polarized, when the angle of the line of 
sight to the magnetic field d is different from the pitch 
angle. Relativistic beaming requires this difference Ad to 
be smaller than 7 . The left- and right-handed CP for 

different signs of Ad nearly cancel for an isotropic particle 
distribution, with a residual proportional to i / -1 / 2 cosd, 
and cancels completely for a pure electron-positron popu¬ 
lation. Any deviation from an isotropic pitch angle distri¬ 
bution will enhance the resulting circular polarization. 
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2.2. Radiative Transfer 


The transfer of polarized synchrotron radiation in ho¬ 
mogeneous astrophysical plasma was derived by Sazonov 
( |l969| ) in terms of nearly transverse electromagnetic waves 
(see Appendix ^). The transfer equations can be formu¬ 
lated for t he four Stokes parameters /, Q , U, V (e.g., Jones 
& O’Dell 1977a| ). The effects included are emission and 
absorption for /, Q. V separately^. The derivation of the 
absorption and general rotation coefficients are given in 
Appendix B| and |^. For power-law distributions of elec¬ 


trons (Eq. D.l): N(E) oc 7 s between 7 m i n and 7 max , and 


< V < 7max 


max? 

ub, the transport coefficients are sum¬ 


marized in Appendix jDj. Here one has to be aware of the 
ultra-relativistic limit used in the derivations. They are 
only useful for y/Jsint? 1. This requirement is not al¬ 
ways fulfilled when the power-law extends down to 7 of 
a few, as will be the conclusion for Sgr A*. In addition 
the shape of the distribution at these energies and the as¬ 
sumed perfect isotropy of the pitch-angle distribution are 
uncertain. The principal treatment of more general dis¬ 
tributions is presented in the Appendix. Linear polarized 
emission is a fixed fraction of the total emission and the 


relative emissivity of CP (Legg & Westfold 1968) is 


Woc(-r 1/2 cot <?, 


Vi 


r'B 


( 2 ) 


where i9 is the angle between the line of sight and the mag¬ 
netic field, which equals the pitch angle for electrons radi¬ 
ating into the direction of the observer. Due to Kirchhoff’s 
law this relation also holds for the absorption coefficients 
Kv/isi = Vv/v /• Normal modes in a magnetized plasma 
are generally elliptically polarized. They are circular for 
propagation along the magnetic field and coupling to the 
gyration of electrons induce different refractive indices for 
left- and right-handed modes. This leads to Faraday ro¬ 
tation described by the transfer coefficient kf, which has 
the dimension of an absorption coefficient, so that after 
propagating a path length A l linear modes rotate by an 
angle = Al kf- 

Due to the steep frequency dependence of all rotation 
coefficients Faraday rotation is important at low frequen¬ 
cies 


kf oc 



-2 


lnyn 


7 s "! 1 

/min 


■ COST? 


( 3 ) 


with the power-law index s for the electrons. Intrinsic 
Faraday rotation is dominated by the low-energy end of 
the electron population and depends on the line of sight 
average parallel magnetic field in the source. 

Perpendicular to the field the normal modes are lin¬ 
early polarized and different refractive indices for modes 
parallel and perpendicular to the S-field leads to bi¬ 
refringence and cyclic transformation of U into V with 

1 A suitable transformation of coordinates leaves U parallel 

to the local B-field. Emission and absorption for U vanishes in 

that situation. 


a rotation angle sin A(/ = A l kq when propagating a 
distance Al. This effect was termed ’repolarization’ by 
Pacholczyk (1973) and is described by the transfer co¬ 
efficient kc (see Appendix 0 
-3 


"B 


kc oc — ) sin 1 2 0'y ] 


( 4 ) 


if the source is self-absorbed and emission at frequency v 
is dominated by high-energy electrons with 7 ra( j 7 min 
or equivalently v vb !^If we have additional cold 
plasma in or surrounding the emission region, the radia¬ 
tion is modified along its path by additional Faraday ro¬ 
tation and conversion k^ and k'q (see Appendix |cj). It 
must be noted that no substantial amount of cold plasma 
ri( c ) « ri( r ) compared to the density of relativistic parti¬ 
cles can exist in linearly polarized, compact synchrotron 
sources ( Tb ~ 10 10 ' 12 K) without depolarizing the emis¬ 
sion (Jones & O’Dell 1977b| ). Inside a source with a homo¬ 
geneous and well ordered magnetic field Faraday rotation 
is necessary for conversion, because no linear polarized 
component U with electric field parallel to the stationary 
B-Field is otherwise produced. 

An electron-positron plasma on the other hand shows 
no Faraday rotation and changing magnetic held direc¬ 
tions along the line of sight either in an ordered held struc¬ 
ture and/or in a turbulent held is required to start conver¬ 
sion. The contribution of protons to Faraday rotation and 
conversion is weaker by a factor (m e /m p ) 3 and ( m e Jm p ) 2 
respectively. In view of the 7 m i n dependence in Eq. (g) this 
implies that Faraday rotation by cold protons is as impor¬ 
tant as by electrons, if 7 min/(I n 7 min) ~ 3/4(TO p /m e ) 3 . So 
for s = 2 a 7 m i n ~ 475 corresponds to equally strong 
rotation by electrons and protons and consequently no ro¬ 
tation. This is not a strong constraint, but for s = 3 we get 
equally strong rotation from cold protons and relativistic 
electrons for 7 m i n = 50 and for larger 7 m j n the protons 
dominate rotation. 


2.3. Intuitive Approach 

To visualise the effects and properties of radiation transfer 
in a compact medium and achieve a more intuitive under¬ 
standing of the relevant effects on polarization, we also 
provide a more qualitative discussion here. For simplicity 
let us separate Faraday rotation from conversion and only 
picture purely linearly or circularly polarized waves in a 
homogeneous magnetic held. 

The two orthogonal normal modes for propagation per¬ 
pendicular to the magnetic held are linearly polarized and 
a purely circularly polarized wave is split into the two nor¬ 
mal modes with a relative phase shift as shown in Fig. 
Without a phase-shift the wave will be purely linearly po¬ 
larized. If, for example, a locally homogeneous magnetic 
held vertically pervades the box in Fig. [i] along the 2 - 
direction, electrons or positrons will be free to move along 
the held lines and resonate with the vertical mode but 
hardly resonate with the horizontal mode along the in¬ 
direction. This yields the bi-refringence discussed above. 
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Fig. 1. A circularly polarized wave can be composed of 
two orthogonal linearly polarized modes shifted in phase. 
A phase shift would be produced by a plasma in a mag¬ 
netic field perpendicular to the propagation direction of 
the waves (here along the 2 -direction). Without phase- 
shift the sum of the two modes would be a purely linearly 
polarized wave. The accompanying movie shows the effect 
of how phase-shifts in a region will turn such a linearly 
polarized wave in to a circularly polarized wave (conver¬ 
sion). 


The resonating electrons or positrons will themselves act 
as antennas and emit a somewhat delayed wave that inter¬ 
feres with the incoming vertical mode, leading to a slight 
phase-shift between vertical and horizontal mode. The ef¬ 
fect of this shift is shown in the accompanying animation^] 
of Fig. [I], where the resulting wave is circularly polarized 
and switches from linear to circular polarization as a func¬ 
tion of the shift. 

Conversion acts also on initially only linearly polar¬ 
ized radiation. The amount of this conversion will depend 
on the misalignment between the incoming wave and the 
magnetic field direction since, obviously, a phase-shift be¬ 
tween two orthogonal modes will have little effect if one 
mode is very small or non-existent. Moreover, a random 
distribution of magnetic field lines on the plane of the sky 
will reduce circular polarization from conversion in exactly 
the same way as linear polarization would be reduced. 

Analogous to the picture for conversion, one can view a 
linearly polarized wave as composed of two circularly po¬ 
larized normal modes when propagating along the mag¬ 
netic field. This is sketched in Fig. where we will as- 


2 See also the authors we bpage 
bonn.mpg.de/staff/hfalcke/CP| 


at 


http: //www.mpifr- 


Fig. 2. A linearly polarized wave can be composed of two 
orthogonal circularly polarized modes shifted in phase. A 
phase shift would be produced by a plasma in a magnetic 
field along the propagation direction of the waves (here 
along the (/-direction). The accompanying movie shows the 
effect of additional phase-shifts on the linear polarization, 
leading to Faraday rotation. 


sume a longitudinal magnetic field, i.e. a field along the 
//-direction. The circular modes will resonate with either 
electrons or positrons gyrating around the magnetic fields. 
The latter will again emit a circularly polarized wave, pro¬ 
ducing a phase-shift when interfering with the incoming 
wave. The effect of the phase-shift in the circular modes 
is shown in the accompanying animation of Fig. where 
one can see that the resulting linearly polarized wave is 
simply (Faraday) rotated. 

An important conclusion to remember therefore is, 
that conversion is mainly produced by magnetic field com¬ 
ponents perpendicular to the line-of-sight or photon direc¬ 
tion, while Faraday rotation is produced by magnetic field 
components along the line-of-sight. Moreover, one can also 
see that conversion is insensitive to the electron/positron 
ratio while Faraday rotation is not. In Fig. [1] an electron 
and an positron are both free to move along the 2 -axis. 
While they will respond in opposite directions to the in¬ 
coming wave, their respective emitted waves will also have 
opposite signs because of opposite charges and hence be 
identical. In the case of Faraday rotation, the incoming 
left- or right-handed circularly polarized wave will only 
resonate with the particle that also has the correct hand¬ 
edness in its gyration - either electron or positron depend¬ 
ing on the magnetic field polarity. A pure pair plasma 
would therefore produce exactly the same phase shift in 
left- and right-handed modes and not produce any net 
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Faraday rotation. In the case of a charge-excess, the di¬ 
rection of Faraday rotation depends on the sign of the 
charge-excess (presumably electrons) and the polarity of 
the magnetic field. This will indirectly also affect the sign 
of circular polarization, if Faraday rotation is the ultimate 
cause of the misalignment between the plane of polariza¬ 
tion and the magnetic field direction. 


3. Outflow/Jet Models 

3.1. Conical Outflows 


Models of flat spectrum radio cores in AG N ass ume in 
general a co nical jet (Blandford & Konigl 1979 ; Falcke 
& Biermann 1995 ), in which plasma is flowing out with 
constant velocity v = j3c and constant half opening angle 
9 = arcsin {R/z), where 2 is the coordinate along the jet 
and R the local radius of the jet. The magnetic field in the 
jet must have an ordered component, which leads to per¬ 
sistent polarization, and probably a turbulent field])]. The 
ordered large-scale field can be separated in modes, which 
carry magnetic flux and will therefore decay as ~ z ~ 2 , and 
modes without magnetic flux (e.g. toroidal fields), which 
behave as 


Otherwise a cosmic conspiracy of electron cooling (adia¬ 
batic and radiative losses) and geometric changes, which 
includes the evolution of magnetic fields, must be pro¬ 
posed. The variety of observed spectral indices between 
1 and 20 GHz indicates, that intermediate stages, where 
several of the effects are present, are quite common. The 
particle distribution without pair-production N(E , z) in a 
cone can vary as 

/ \ -2(l+a/3) 

N(E,z)=N(E,z in )(^j (6) 

where N(E, Zi„) = N(E) is the energy distribution at the 
injection point. The parameter a is zero for a freely ex¬ 
panding and over-pressured jet and a = 1 for adiabatic 
losses in pressure equilibrium. Radiative cooling can fur¬ 
ther change the distribution and may be considered sepa¬ 
rately. For a conical outflow with partial adiabatic particle 
cooling the resulting sychrotron spectrum can be sepa¬ 
rated in three spectral regimes. 

{ a = 5/2 optically thick 

ol = | (l - self-absorbed (7) 

a=(s — l)/2 optically thin 


B = B 0 



(5) 


The second type of modes dominate at large radii. One 
possible configuration is a force-free equilibrium for a par¬ 
allel jet (Konigl & Choudhuri 1985| ). The resulting struc¬ 
ture can show both a mode with magnetic flux and modes 
of twisted flux tubes with toroidal and two oppositely di¬ 
rected poloidal fields, so that the magnetic flux along z 
vanishes for the second modes. Other configuration may 
not be force-free but are equally valid. The most likely 
structure is a helical R-field due to the rotation of the 
footpoints in an accretion disk. The spiral can be re¬ 
garded as the superposition of a field component B z along 
z and a toroidal component Btf,. The slope of the spiral is 
as = arcsin (B^/Bz). 

Every jet in perpendicular pressure equilibrium (re¬ 
direction) with its surrounding will suffer adiabatic en¬ 
ergy losses. Relativistic electrons, which are injected at the 
base of the jet, will cool down due to adiabatic expansion, 
which leads to inverted radio spectra observed in some 
core dominated extragalactic radio sources. Further along 
the jet electrons have to be reaccelerated. Alternatively 
the jet is highly over-pressured relative to its surround¬ 
ing and adiabatic losses are negligible consistent with flat 
radio spectra for conical jets cores. It should be noticed, 
that toroidal fields will induce electric currents and hook 
stresses will confine the jet, if the magnetic field is well 
ordered and strong enough to influence the jet dynamics. 


3 Alfvenic turbulence is expected from diffusive particle ac¬ 
celeration at shock fronts. Pitch angle scattering at Alfven 
waves will also lead to isotropic particle distributions, which 
we assume here. 


The optically thick regime can only be observed if the out¬ 
flow terminates at a maximal distance, or fragments into 
subcomponents, which break the self-similarity of the con¬ 
ical model. In the self-absorbed regime, the flux is domi¬ 
nated by emission from the region around the optical sur¬ 
face (r=l-surface). A sufficiently large ratio z ou t/z- m of 
outer to inner radius of the self-similar conical outflow is 
necessary for it to be observable. Optically thin emission 
is always present at frequencies above = ^abs^in), if 
not other radiation processes like bremsstrahlung or dust 
emission dominate. In some sources the low-frequency 
spectrum can change due to free-free-absorption. 


3.2. Turbulence in Jets 


Synchrotron emission from self-absorbed radio sources 
with brightness temperatures Tb of ~ 10 n K imply 
near energy equipartition of radiating electr ons and/or 
positrons and magnetic field (Readhead 1994 ). This find¬ 
ing is modified by relativistic Doppler boosting for vari¬ 
able flat spectrum radio cores. VSOP-observations show 
that the observed Tg can be larger than 10 12 K (Bower & 
Backer 1998; Tingay et al. 2001) in selected sources. This 
agrees with the observed superluminal speeds /3 0 bs ~ 5, 
which require Doppler factors > 5. The fractional LP of 
radio cores is usually only a few per cent or smaller and re¬ 
quires either strong Faraday depolarization (Tribble 1991) 
or tangled B— fields. The observed rotation measures in 
some quasars (Taylor [2000 ) are not sufficient to depolar¬ 
ize the radio emission at cm- wavelength]^] and we infer the 
presence of a turbulent contribution to the global R-field 


This is not true anymore, if the rotation measure rises 
sharply towards the centre {RM oc 2 -2 or steeper). 
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described above. The turbulent field is effectively depo¬ 
larizing the source, if the amplitude of the dominating 
Fourier modes of the turbulent field are about a factor of 
2 larger than the local contribution of the global field, so 
that local field reversals occur. We describe the turbulent 
magnetic field as a local superposition of incoherent waves 
with wavenumber k, which decorrelate over distances of a 
few times their wavelength. For the treatment of radiative 
transfer in turbulent jets, we consider the turbulence to 
be frozen in and time evolution to be unimportant. 

The turbulent wave spectrum is characterized by an 
outer wavelength and corresponding wavenumber fc out = 
27 t/L and a dissipation wave number kd ■ Between these 
wave numbers an inertial range with energy cascading 
from small to large wave numbers will develop^]. The local 
strength and orientation of the turbulent magnetic field 
will be determined by modes with wave numbers around 
fcout and the typical length-scale for changes of the mag¬ 
netic field is \B/VB\ ~ k~^ t . 

We assume that the spectral energy density in the in¬ 
ertial range is described by a Kolmogorov spectrum 

E(k) = F e 2/3 fc" 5/3 . (8) 


Here F is a number of order unity. The dissipation wave 
number is related to t he en ergy dissipation rate e by the 
inequality (e.g. Frisch 1995 ) 


hi > e v ~ 3 


(9) 


where v is the viscosity. The dissipation e for Alfvenic 
turbulence is then given by 


= ( B 2 T /(B 2 0 F)f 2 k out v\ 


v 2 a = 




-i£- (10) 


with va the Alfven velocity in the global field and Bj, / B q 
the energy density ratio of turbulent to the global field. We 
can combine the dissipation rate © with the estimate for 
the dissipation wave number (|j) to get the viscosity 

v > (B^/B^F^ik^t/kd) 1 / 3 VAk- 1 . ( 11 ) 


If the energy dissipation process is known, it is possible to 
estimate the thermalisation of magnetic energy along the 
jet. One possibility is a first resonance with the gyration 
of protons or relativistic electrons. 


3.3. Error Estimate for Radiative Transfer in Turbulent 
Fields 

A numerical treatment of radiative transfer in turbulent 
magnetic field will not resolve incoherent fluctuations with 
wave numbers larger than k n depending on the numerical 
code. The largest effect possible of the unresolved modes 
with k n < k < kd is a coherent addition. The critical wave 

5 Here it is implied, that an instability exists, which injects 
energy into the turbulent cascade at the outer wavelength. 



Fig. 3. Outflow model for the radio spectrum of Sgr A*. 
The result of model calculations for total flux I (solid line), 
linear (dense shaded area) and circular polarized flux V 
(sparse area) are shown for a distance of 8 kpc. The shaded 
areas mark the expected variability due to turbulence with 
fcout-R = 50. The global magnetic field structure is a spira l 
with B^/Bz = 1. CP data are taken from Bower ( [2000 ), 
averaged from two campaigns. LP measurements are from 
Bower et al. ( |L999c ) and (2001) 


number fc c , for which the rest k > k c can result in field 
reversals is 


kn = kn 


bout B ko 


1 


B 0 6 < cos cj) k > k 


+ 1 


( 12 ) 


We expect to overestimate emission and absorption of po¬ 
larized radiation, which depends on the projection of the 
magnetic field onto the line of sight and which change sign 
with field reversals. Using the estimate of < cos 4> k > k ~ 
2/7r, we have an upper limit to the reduction of the effec¬ 
tive absorption and emission coefficients by 

2 6B 0 

p = -~- 

AH C 7r fcout B k o 



In the numerical models the smallest length-scale k n 1 is 
adjusted, so that p > 1 and no reduction has to be applied. 


4. Depolarization and Conversion 


4.1. Analytical Estimates 


Turbulence leads to a reduction of measured polarized 
flux, if observations are not able to resolve the largest 
turbulent scales L m kffut i n the flow. Additional reduc¬ 
tion of polarization occurs along the ray path through the 
source, if the Faraday optical depth in regions of size L is 
large. We define the Faraday cell depth as 


Tp — rtp/kout : 


(14) 


where up is the Faraday rotation coefficient given in (D._3). 
At a given frequency only the polarized flux from the opti¬ 
cal surface is important and will be considered, even if the 
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Fig. 4. Fractional polarization for the numerical calcula¬ 
tions ( 7 tl- boxes; ire'- crosses) shown in Fig.|] compared 
to the upper limits for linear (arrows) and the measured 
circular polarization (triangles). Data are taken from two 
campaigns (Bower 2QOO| ). The Percentage of linear polar¬ 
ized flux (solid line) and circular (dash-dotted) for the 
analytic estimates of Eq. (|l9|) and Eq. ©) is also shown. 


fractional polarization in the optically thin region above 
the optical surface is larger. Both conversion and rotation 
changes the linear polarized flux and their relative impor¬ 
tance is measured by 


£ 


KF 

KG 


s + 2 v cos "d In 7 m i n 
S+Ivb sin 2 1? Train 1 



.(15) 


For flat spectrum self-absorbed outflows the ratio v/vb = 
T'rad is nearly constant and radiation is dominated by elec¬ 
trons with the same 7 ra d at all radii. For the special case 
s = 2 we get 


^ _ 4 cos 1 } / 7 rad \ 2 In 7 min 

3 sin $ \ Tmin ) Tmin lll(7rad/ Tmin) 

The dependence of £ on the viewing angle 1 ? is £ oc 
cos $/ sin 2 d and the ratio £ is shown in Fig. || as a function 
of brightness temperature Tb oc Trad for different values of 
7min ■ Faraday rotation is less important than conversion, 
if Trad > 20 and Tmin ~ 0.97 ra d- Whenever the power-law 
population of electrons extends below 7 m i n = 20, Faraday 
rotation is always stronger than conversion with the ex¬ 
ception of almost perpendicular magnetic fields. For high 
brightness temperatures and 7 m i n < 0.57 ra d Faraday ro¬ 
tation depolarizes the emission, which is applicable to the 
situation in Sgr A* discussed below. 

The polarization is determined close to the optical sur¬ 
face and the relative Faraday optical depth is 


2(s + 2) / TVad V In Tmin 

(^ T l)^out-^ \Tmin/ Tmin 


(17) 


Fig. 5. The ratio £ of Faraday rotation to conversion as 
function of brightness temperature for the intermediate 
angle $ 0 , defined by cos$o = sin 2 $ 0 - The electron spec¬ 
tral index is s = 2.5 and £ is plotted for four values of 
Tmin =2 (dashed line), 5 (solid line), 10 (dash-dotted), 
and 20 (dashed-3xdotted). 


where k ou tR is the number of turbulent cells along the ray 
path. In analogy the relative strength of conversion is 

Tc/r ~ (.s - 2) fcouti? 


Trad 

Tmin 


s -2 


- 1 


(18) 


Consequently for Trad ~ 100 substantial Faraday rotation 
occurs for 7 min less than 10. 

In sources where the relative Faraday optical depth 
© is larger than unity the depolarization along a ray 
path is dominated by Faraday rotation within cells of 
size fc c 7 ut . The optical surface of the source is covered by 
(k out R) 2 cells. The instantaneous fractional LP due to in¬ 
ternal Faraday rotation is then 


ttl = 


a/Q 2 + u 2 s + 1 1 


s + 7/3 fcout R 


t/t f . 


(19) 


The time average ( 7 r^) will vanish in a stochastic magnetic 
field. In jets or outflows from rotating central objects we 
expect an ordered magnetic field component. For the spiral 
magnetic field described in Sec.3.1 the averaged field in 
the plane of the sky is B z sin 7 ? and we expect a mean 
fractional LP of 


(tt) 


s + 1 


1 


By 


s + 7/3 koutR t f B 0 


sin#. 


( 20 ) 


The field strength Bq is the sum of globally ordered and 
turbulent field and 6 the angle between the jet direction 
and the line of sight. 

The appearance of linear polarized radiation perpen¬ 
dicular to the local magnetic field is the starting point 
for conversion to circular polarized radiation. Linear po¬ 
larized radiation from outside the cell will have a non¬ 
vanishing Stokes U locally, while LP emission intrinsic to 
the cell must undergo Faraday rotation before conversion 
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Fig. 6. The same model for Sgr A* as shown in Fig. 1 
with a tightly wound spiral structure B ( f > /B z = 3 seen 
at an angle of 30° to the jet axis. It can not explain the 
observed CP. 


can take place. In the case of dominating Faraday rota¬ 
tion Tp > t the fraction of suitable LP along a particular 
line of sight is (p + 1 )/(p + 7/3)t/t f . This fraction can 
be converted with an efficiency tc/t , if tq < r. The re¬ 
sulting fractional CP follows from an average over surface 
elements of relative size l/(fc 0 ut-R) 2 : 


vc = ^ (j—n) ( t/t f )(t c /t ). ( 21 ) 

I s + 7/3 \k out RJ 

Again the time average CP in a stochastic field will van¬ 
ish ( 7 re) = 0. In an outflow with a spiral magnetic field 
component the average CP is approximately 


(tt c ) 


0.5 


s + 1 
s + 7/3 



tc B z 

-— cos 6 . 

tf Bo 


( 22 ) 


If Faraday rotation within one cell is small, the rele¬ 
vant perpendicular polarized emission is proportional to 
(fcout-ft ) -1 / 2 an d the resulting CP 


7TC 


1 


/•’out R 


3/2 


(tc/t) 


(23) 


In that case LP is not reduced by Faraday depolarization 
along a particular line of sight and 

(24) 

p -\- 7/3 y rC ou t R ) 

the resulting linear polarized flux is larger than the circu¬ 
lar by a factor t/tc, which is larger than one. 


4.2. The Specific Model for Sgr A* 

The inverted spectrum of Sgr A* implies either a non- 
conical geometry of the outflow, a magnetic field com¬ 
ponent, which decays faster than 1/z along the outflow, 
or reduced, but not absent, adiabatic cooling of the radi¬ 
ating, relativistic particles. Acceleration of the ouflowing 


plasma (Falcke 1996 ) also produces inverted radio spec¬ 
tra. For simplicity, this is not explicitly considered here, 
however, its effect is essentially covered by the assumed 
scaling of magnetic field and density. Cooling of relativis¬ 
tic particles implies IV( 7 , z) oc 7 ~ s z ~ 2 i 1 + a /3) ^ where a = Q 
for freely expanding outflows, and a = 1 for adiabatic 
losses in pressure equilibrium with the surrounding gas. 
The spectral index a = 0.3 implies a = 0.1(s + 4) ss 0.7 
for our choice of the electron spectrum s = 3. This is one 
of the possible parametrisations of the observed inverted 
spectrum. The magnetic field at the base of the outflow is 
60 G at Zi n = 10 Rs in our model, where the width of the 
outflow equals the distance from the black hole of mass 
2.610 6 Mq. The presumed electron distribution extends 
from 7 min = 5 to 7 max = 250. The high energy cut-off is 
required by the infrared limits of Sgr A* spectrum and re¬ 
sembles a wide quasi-monoenergetic distribution (Beckert 
& Duschl 1997). The low-energy end is set by Faraday 
rotation and conversion to produce the observed polariza¬ 
tion in an inhomogeneous, optically thick jet or outflow. 
The emission becomes optically thin above v = 5-10 11 Hz, 
which implies that optical depth is unity for 7 ra( j « 35 at 
the base of the outflow. Due to the inverted spectrum 7 ra( j 
varies with radius 7 ra d oc --° 068 j which has to be consid¬ 
ered in Eq. 0 and Eq. H This implies a weak frequency 
dependence of (ttl) and (no) of the model seen in Fig |[ 
Because Faraday rotation is strong and conversion within 
one cell is weak in our model with th/t ~ 6.04/(fc out I?), 
where we use an outer scale k ou t = 50/1?, we can use Eq. 
( p2[ ) and Eq. ( f20| ) to estimate the mean CP and LP. The 
analytical estimates for ( 717 ) and (ire) are shown in Fig. 
[| together with the measured CP and the results from 
numerical solutions of the transfer problem. These esti¬ 
mates hold as long as the outflow is self-absorbed below 
3 • 10 n Hz. 


The model for Sgr A* requires depolarization domi¬ 
nated by a large Faraday rotation depth. In this case the 
outer turbulent scale fc ou t is poorly constraint. For the 
analytic model of Eq. (j|(|) and (|||) to be valid fc ou t has 
to satisfy 6 <C koutR 1250. The lower bound comes 
from tc/t -C 1 and the upper bound from t f /t 3> 1 . 
The value for the numerical treatment shown in Fig. and 
Fig. [| is k out R = 50. The resulting spectra for numerica l 
solutions of the radiative transfer problem (Eq. A.5-A.8) 
on a 100 2 grid covering the jet seen under an angle of 
9 = 30° in the rest frame of the gas is shown in Fig. 
The required electron density is 2.8 10 7 (z/10 Rs)~ 2 cm " 5 
in a global H-field of 60 (z /10 i ?^)^ 1 G. The half opening 
angle of the subsonic outflow is 4.5° with a bulk motion 
of /3 = 0.4. It turns out that a spiral structure for the 
global magnetic held seen under an angle 9 < 7t/2 — as is 
preferred for reproducing the level of linear and circular 
polarization. The numerical simulation use a as = 7 t/ 4 - 
spiral corresponding to B z /Bq = 1 / x/ 2 . In the limit of a 
very long spiral with as 0 the jet must be seen at an¬ 
gle smaller than ~ 40°, because conversion and Faraday 
depolarization is stronger for small 9 in our model. For 
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Fig. 7. The spectrum of a typical inner jet component 
in total intensity I (solid line), LP (triangles), CP (dia¬ 
monds). The fluxes are normalized to the peak flux in I. 
Energy equipartition between electrons and 13-field is as¬ 
sumed. The electron spectral index is s = 2.5 and the 
power-law extends from 7 m ; n = 20 to 7 max = 10 4 . The 
component is not homogeneous but is part of a jet with 
Zout/zin = 10 and a tight helical field ( B z /B$ = 1/3) seen 
under an angle of 55° in the rest frame of the radiating 
plasma. 


the assumed electron-proton plasma the kinetic bulk en¬ 
ergy is about 7.4 1 0 39 erg/s. Together with the magnetic 
energy flux, the thermal energy flux derived from the 
half opening angle of the flow, and the energy supply 
needed to overcome the gravitational potential with the 
large mass loading of the flow starting at 10 Rs, the total 
power is Lj e t ~ 1.3 10 4 ° erg/s. From Eq. (0) we can es¬ 
timate the turbulent energy dissipation rate along the jet 
to be 2.7 1 0 39 ln(z out /2in)[(fcout^)/50] erg/s. The inverted 
part of the spectrum of Sgr A* extends from 1 GHz to 
350 GHz and the ratio z ou t/zi n must therefore be larger 
350. This provides a upper bound to the outer turbulent 
scale k out R < 24 so that the kinetic energy is not dissi¬ 
pated before reaching z oui . The total jet power is 5 or¬ 
ders of magnitude larger than the emitted radio luminos¬ 
ity, which increases the required accretion rate to power 
the jet. The brightness temperature of the model is little 
less than the equipatitition temperature and the inverse 
Compton luminosity from optical to X-rays is much less 
than the radio luminosity. The radio jet of Sgr A* is a very 
inefficient radiation source. For the polarization in M81* 
we would qualitatively obtain similar numbers, with how¬ 
ever, a higher jet power. 


4.3. Jet Components in Quasars 

In the brigh t jet sources 3C 84 and 3C 273 (Homan & 
Wardle 1999 ) it has been demonstrated, that CP can be 
detected in the core and the innermost jet component. The 
degree of LP is equal or less than CP in the inverted spec¬ 


trum cores. Various other components show 0.5% circular 
polarization in these sources. For 3C 273 it is claimed, that 
the circular polarization is predominately associated with 
new ejected jet components. This has to be taken with 
caution as Taylor ( 1998| ) reports a large rotation measure 
of RM = —1900 radm -2 for the core of 3C 273 and almost 
equally large RM for the core of 3C 279. The reduced lin¬ 
ear polarization in 3C 273 may therefore be due to depo¬ 
larization in surrounding cold gas. They also report that 
the component CW in 3C279 is +1.2% circularly polar¬ 
ized, while the linear polarization of component CW is 
13% on average. A model for a typical jet component is 
shown in Fig.y. In contrast to Sgr A* we have here less re¬ 
duction of polarization due to turbulence (as we only look 
at one spatially resolved jet component) and less Faraday 
depolarization. The model invokes a higher 7 m i n of order 
20 and a power law up to 7 ma x = 10 4 that produces emis¬ 
sion well into the optically thin regime, which is known 
not to exist in Sgr A*. This gives one the characteristic 
LP-to-CP ratio of ~ 10 observed in quas ars and recently 
also in X-ray binaries (Fender et al. 2002). 

It is also interesting to note that we can produce 
the observed CP with a 7 m i n of order 20. This is some¬ 
what higher than the rather low values found by Wardle 
et al. ( 199^ ) and does not place quite so stringent con¬ 
straints on the energy budget and the matter content of 
the jet. Somewhat more realistic energy distributions at 
low-energies other than a sharply cut-off power-law may 
further relief these constraints. 


5. Polarization Variability 


In the presence of a turbulent magnetic field the degree 
of polarization (both circular and linear) and position an¬ 
gle will vary stochastically with a timescale At = A z/va- 
Polarization variability is expected to be faster than varia¬ 
tions in total flux, because the relevant length-scale is the 
outer turbulent length-scale k~ ut . With va ~ cs ~ c/V3 
we get a characteristic variability time 

At = x oc j/—( s + 4 )/( s + 4 + 4a /3) (25) 

c(fcout-R) 


The effective radius at a given frequency can be deter¬ 
mined from a model for the total flux, while the turbu¬ 
lent scale can be derived from polarization of optically 
thin emission or estimated from the ratio (ttl) / (ttc)- 
Turbulence also leads to decorrelation of polarization and 
position angle across frequency bands. Polarization prop¬ 
erties from optical surfaces a distance fc“ u 1 t apart should 
be uncorrelated. This translates into a relative distance 
A z/z = sin 6/(k ou tR) and both LP and CP should vary 
across frequency bands An set by the turbulence in the 
source 


Au ( a/3 \ sinf? 

v V + s/4+lJ fc out i? 


(26) 


The suggested model for Sgr A* with the Alfven veloc¬ 
ity of va = 210 9 cm/s and fc ou tR = 24 implies a vari- 
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ability timescale of At sa 16 h at 1 GHz and an accord¬ 
ingly shorter variability time scale at higher frequencies. 
The decorrelation across frequency bands is expected at a 
width Av/v « 2.4 • 10 -2 . Longer integration times than 
set by Eq. (^5|) and frequency bands wider than given by 
Eq. @ will produce smaller polarization measurements 
than intrinsically available in the source. 


6. Discussion 


Recent observations of radio circular polarization in AGN, 
X-ray binaries, and the Galactic Centre black hole suggest 
that CP at the 0.3%-1%-level is common to many self- 
absorbed synchrotron sources. Faraday rotation and con¬ 
version in a magnetized and therefore bi-refringent plasma 
produce enhanced circular and reduced linear polariza¬ 
tion. Both processes are sensitive to the presence of low- 
energy electrons and the orientation of the global magnetic 
field. 

The standard jet model for compact radio cores with a 
helical plus a turbulent magnetic field can well reproduce 
the circular and linear polarization spectrum of sources 
like Sgr A* and M81* with their high CP-to-LP ratio. 
The suppression of LP is achieved by the presence of a 
significant number of low-energy electrons in the source 
and an absence of an optically thin power-law extending 
to higher energies. The same model can also explain the 
typical level of circular polarization in blazars and the 
CP-to-LP ratio observed in blazars and X-ray binary jets. 
In this case a number of low-energy electrons is reduced 
with respect to the Sgr A* model and a power-law in the 
electron distribution exists. 


For Sgr A* the number of low-energy electrons produc¬ 
ing conversion and depolarization needs to be significantly 
higher (by 2-3 orders of magnitude) than the number of 
radiating hot electrons. This means that a large fraction 
of the outflowing jet material is in the form of hidden 
matter shielded by self-absorption. This increases the es¬ 
timates of the total jet power, which can be 5 orders of 
magnitude higher than the synchrotron luminosity. If one 
presumes that this power has to be provided by an accre¬ 
tion flow, the minimum accretion rates of lO -9 ” -8 M 0 /yr, 
previously estimated from “maximally-efficient” jet mod¬ 
els for Sgr A* (Falcke et al. 1993; Falcke & Biermann 


1999) need to be raised to about 10 6 M 0 /yr. This is quite 


consistent with recent estimates of Bondi-Hoyle accretion 
rates onto Sgr A* (Baganoff et al. 2002| ) and with sug¬ 
gestions for a coupled jet plus ADAF model (Yuan et al. 
2002 ), where the emission from the accretion disk is highly 
suppressed with respect to the jet. 

It is also interesting to note that to fit the CP with con¬ 
version one requires an asymmetry in the magnetic field 
components. This is naturally achieved by a helical mag¬ 
netic field as is presumed to exist in jets. A symmetric 
configuration, e.g. a tightly wound helix or even a toroidal 
magnetic field structure would have difficulties to produce 
the observed level of CP. 


The stable handedness of CP over 20 years also im¬ 
plies a long-term stable component of the unidirectional 
field along the line-of-sight. This indicates that the polar¬ 
ity of the magnetic field (the “magnetic north pole”) has 
remained constant over the last two decades. In view of the 
rather short accretion time scale in Sgr A* one could also 
speculate that this polarity is related to the accretion of a 
stable large-scale magnetic field which is accreted and ex¬ 
pelled via the jet. The same can be said about blazars and 
X-ray binaries, where the stability found in GRS1915+105 
by Fender et al. ( [2002 ) is particularly interesting since the 
intrinsic accretion time scales in X-ray binaries are much 
shorter than those in supermassive black holes. 


Another important aspect of CP measurements is the 
question of the matter content of jets. We find that the 
constraints from CP of individual jet components for the 
jet power in blazars are not quite as severe as previ¬ 
ously claimed and a statement in support of a pure elec¬ 
tron/positron jet has to viewed with caution. For Sgr A* 
or M81* the situation may be different. If the depolar¬ 
ization is indeed i ntrins ic to the jet and not a surround¬ 
ing medium (Agol 2000| . Quataert & Gruzinov 2000 ), one 
needs a high Faraday optical depth in the source, which 
can only be achieved by an excess of “warm” ( 1 ^ 7 ^ 
100) electrons in an electron/proton plasma. 


While we have here assumed that all electrons are dis¬ 
tributed in a single power-law, the actual situation may 
be quite different. For Sgr A* a power-law is actually not 
needed and we could obtain rather similar results with a 
two-temperature electron distribution, with temperatures 
corresponding to 7 m i n and 7 max respectively. This is not 
quite possible in blazars or bright X-ray binary jets, where 
extended electron power-laws are directly observed in the 
optically thin regime. It could well be that the radiative in¬ 
efficiency of Sgr A* is due to the lack of effective shock ac¬ 
celeration that would increase the number of high-energy 
electrons with respect to the number of low-energy elec¬ 
trons (and in turn decrease the CP-to-LP ratio). The ori¬ 
gin of these different electron distributions and their role 
for the radio-loudness of jet sources should be a very ex¬ 
citing question for further research. 


By improving our sensitivity and imaging all four 
Stokes parameters at multiple frequencies in the future, 
it will be possible to construct models of the entire emis¬ 
sion and transfer processes in the source and determine 
the composition and energy spectrum of the relativistic 
plasma within jets. 


Appendix A: The Radiative Transfer Problem 

In a weakly anisotropic medium the propagation of elec¬ 
tromagnetic radiation is determined by the dielectric ten¬ 
sor eij . The normal modes in a magneto-active, anisotropic 
plasma are quasi-transverse, but they are not orthogonal. 
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The transfer of the intensity tensor Jy along a ray pat 
k/k is given by (Sazonov 1969, Zheleznyakov et al. 197' 
the emissivity Sij and e^ 


d L 


1] __ Q 

ds ~ bii 




(A.l) 


For transverse waves the intensity tensor is a 2 x 2 tensor 
perpendicular to the wave vector k and for the emissivities 
and the dielectric tensor only the transverse components 
enter Eq. (A.l). The f indicates the complex conjugate 
and transpose of e. 

In a magneto-active plasma with the local magnetic 
field B unperturbed by the wave along the z-axis the di¬ 
electric tensor has only one symmetry e xy = —c yx - The 
tensor e may be separated into a hermitian ( H ) and anti- 
hermitian (A) part according to 


e ij 

= e (H) +e {A) 

13 % 3 

(A.2) 

e (H) 

13 

= 2 ( e b + e ij) 

(A.3) 

e {A) 

13 

= 2 ( e b — e ij) 

(A.4) 


dJ 
ds 
d Q 
ds 
dU_ 
ds 
dV_ 
ds 


= rji ~ K il ~ k qQ - K vV ~ kuU 
= VQ - k iQ - k qI - k fU - HqV 
= T]u — KlU ~ Kjjl + KfQ — KcV 
= T] V ~ KiV ~ Kyi + HqQ + KC'U 


(A.5) 

(A.6) 

(A.7) 

(A.8) 


Ku = - —sin 0(Q(e yz ) + 5(e zy )) . 


(A.12) 


6 In this section |X] the path length along the ray is called s 
not be confused with the spectral index of a power-law distri¬ 
bution. 


The rotation coefficients for Faraday rotation kf and con¬ 
version Kc, h,Q are 

kf = y c (2cos6'3(e xy ) + sin0 (S(e yz ) - S(e ZJ/ ))} (A.13) 


kc = ~y c {^( e w) “ cos 2 05ft(e xx ) - sin 2 09t(e zz ) 
- sin 0 cos 6 (3?(e zx ) + Sft(e xz ))} 


hQ = — sin# (SR(e yz ) + 3f(e ZJ/ )). 


(A. 14) 
(A.15) 


The hermitian eS H ' 1 descri bes absorption processes, while 
the action of in Eq. (AT) conserves the total inten¬ 
sity and rotates the polarization vector of elliptically po¬ 
larized radiation on the Poincare sphere, which is formed 
by the nor malize d Stokes parameters (Q, U,V)/I (Kennet 
& Melrose 1998). This generalized rotation consists of 
Faraday rotation and ordinary conversion between U and 
V and extraordinary conversion between Q and V. From 
Eq. (A.l) the transfer of the Stokes parameters for polar¬ 
ized waves follow 


In the case of an isotropic distribution of unperturbed 
particles f(p) two additional symmetries of the dielectric 
tensor e appear e yz = —e zy and e xz = e zx . Therefore the 
emissivity and absorption coefficient for U and the ex¬ 
traordinary conversion hQ vanish. 

Appendix B: Dielectric Tensor and Plasma 
Kinetic Theory 

The reaction of a distribution of charged particles fo(p) 
in a magnetic field B to a perturbing wave can be 
derived from the linearized first-order pertur batio n to 
the Vlasov equation (Montgomery & Tidman 1964 ). In 
Fourier-Laplace space the perturbation of the particle dis¬ 
tribution fi is given in terms of a propagator G(</') 


fi = ^ d(j)' G(</>') yE — i— A (k A E) 

H J± oo V S 


dfojp) 

dp' 


.(B.l) 


Here E is the perturbed electric held of the wave, e the 
charge and v' the velocity of particles characterized by 
a phase < j>' along their path, s = so — iu> is a complex 
frequency and fi = Ho/7 the gyro-frequency of particles 
with energy qmc 2 . The sign of the lower boundary in Eq. 
(B.l) is determined by the charge of the particles (—oo for 
electrons; +00 for positrons). The propagator is 


G{(j)') = exp 


[ d cf>" (s + ik ■ v")/Cl 

J d>' 


(B.2) 


The transport is described by the transport coefficients 
for absorption 

ki = y c {^( £ to) + cos 2 6 ^s(e xx ) + sin 2 03(e zz ) 

-sin0 cos0(3(e zx ) +3(e xz ))} (A.9) 

kq = ^ {-3(e TO )+ cos 2 0 3(e xa; ) + sin 2 0S(e zz ) 

— sin 9 cos 9 (3(e zx ) + 3(e xz ))} (A.10) 

k v = (2cos 6^R(e xy ) + sin 0 ($t(e yz ) - K(e zy ))}(A.ll) 


which can be used to calculate the Fourier-Laplace trans¬ 
form of the current density in response to the perturbing 
wave 


j = en 0 / d 3 pfi(p)v. 


(B.3) 


where ng is the particle density in the plasma. The com¬ 
ponents of the dielectric tensor can then be read off the 
relation 


djEj = 47 t iuj 


ji , (B-4) 

For isotropic d istrib utions fg ( p ) the second term in paren¬ 


theses in Eq. (BT) vanishes and the calculation of cq is 


greatly simplified. The ^/-integration in Eq. (BT) results 
in a factor —ifi/(nfi + k\\V\\ — to — ie), where e is small 
and positive, because no emission process is described by 
the dielectric tensor. In the limit e —> 0 + the momentum 
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integrals in Eq. (B.3) enclose singularities and must be 
interpreted as 


= V dx 


/(a 


x — z 


kf = 47TCOS0 


-£ 

n= 1 
2 ^02 


cilo 

p'fnil 


V / dp|| / d p± 


nil 


dfojp) 

dp± 


Sc) 


u>pila 

CJjJ 2 


COS0 


(s + 2) ln7„ 


■ cos 9 


Sc) 


t\j /~i — 


' I 2 0 2 

■C p- '() 

2 cu; 3 


sin 2 0 . 




's- 2 


All the absorption and rotation coefficients can be scaled 
to the inverse length-scale 


nif(z ), (B.5) 


lim [ dx /w 

?^o+ J x - [z + ie) 

where V f indicates a principal value integral. The terms 
containing 5-functions lead to the imaginary part of the 
diagonal elements of and to the real part of the off- 
diagonal elements. They describe absorption of the four 
Stokes parameters. The terms connected with principal 
value integrals are 90° off in the complex eij-plane and 
describe the generalized rotation. 

Appendix C: Faraday Rotation and Conversion 


,_i V P v B vb 

L = K 0 = TT -T = T e crl e — 7T , 

c v A 


(D.2) 


with r e the classical electron radius, n e the particle den¬ 
sity, v B = eB / (27ttoc) the cyclotron frequency and v p = 
i/n e e 2 /(7TTO e ) the plasma frequency. The transport coef¬ 
ficient for Faraday rotation by power-law electrons (Eq. 

ED is 


kf = 2 ko cos 6 


s + 2 In 7„ 


S + l 

I min 


(D.3) 


With the help of Eq. (B.5) in Eq. (B.l) we get for Faraday 
rotation (A.12) by an isotropic distribution of particles 

2 OO n p 


and for conversion 


vb . 2 a 

Kc = o—sm 0— 
v s 


2 1 

1 , 

( UB \ 

- 2 1 

U^ 2) 1 



.(D.4) 


” n 2 il 2 - (fc,| V || - w) 2 k±v± Jn(z) 

In general Faraday rotation is represented by an series 
of Bessel-functions J n (z) of order n and their derivatives. 
The argument of the Bessel-functions is z = k±v±/il. In 
the high-frequency cold plasma limit only the n = 1 term 
is important and the momentum integral can be solved by 
partial integration to get the classical limit 


The three absorption coefficients are 


Ki = sm t 


x? ^ r (, n\ r (, V 


4 6/ V 4 6 


kq = 


s + l 
s +7/3 


ki . 


(D.5) 

(D.6) 


(C.2) 


In the ultra-relativistic limit 7 1 Sazonov (1969) de¬ 

rived for power-law distributions in energy with s > — 1 
the Faraday rotation coefficient 

,2f 


a ( UB± \ 

Ky = Ko cos v^- J 


V B ± \ ( s+1 )/ 2 (s + 3)(s + 2) 


av- (s + l) 7 S+1 (C ' 3) 

where the plasma frequency top must be taken for the den¬ 
sity of particles under consideration. In the high-frequency 
limit for a cold plasma the conversion coefficient kc is 


s + l) 
s 7 


3^ u 

4 1412/ V 4 12 


(D.7) 


The emission coefficients can be scaled to the emissivity 

,,2 


?7o = k-S VB m e c 2 = (r e /c) n e u B m e c 2 


(D.8) 


(C.4) 


and for the power-law distribution of Eq. (D.l) we get 


Again Sazonov ( 1969 ) gave an ultra-relativistic approxi¬ 
mation for conversion by a power-law distribution of par¬ 
ticles 

K C(c) _ 0 ( 7rnin ( n n ) 1 •(^■^) 


VI = Vo sin 9 


Vb ± \( s - 1)/2 


3 s / 2 f s 19 

X -;- -F 


,m \ flo sin 9 J 

Appendix D: Transport Coefficients 

Here we summarise the transport coefficients for a rel¬ 
ativistic plasma with a normalized power-law distribu¬ 
tion /V(7)d7 = n e /o(7)d7 above a lower cut-off energy 


.. , r(- - - 

2 (s + l) V 4 12 ) V 4 12 


s + 1 n 
vq = , nl „ vi ; vu = o, 


'Tmin^'e^ 

foil) = ~t~ 77 _S ■ 
7min 


(D.l) 


s + 7/3 

Vv = -Vo COS 9 [— — J 

3(s-i)/2 ( s + 2 ) 

x -^ i 

2s 


(D.9) 

(D.10) 


72 


r, i + f) r u + 5 
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